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ABSTRACT 

Aims. We investigate temperature fluctuations in Hn regions in terms of a two-phase model, which assumes that the nebular gas 
consists of a hot and a cold phase. 

Methods. We derive general formulae for r([Om]), the [Om] forbidden line temperature, and T(Hi), the hydrogen B aimer jump 
temperature, in terms of the temperatures of the hot and cold phases, T h and T c . 

Results. For large temperature differences, the values of t 2 required to account for the observed difference between 7\[Oiii]) and 
T(Hi) are much lower than those deduced using the classical formulae that assume random and small amplitude temperature fluctu- 
ations. One should therefore be cautious when using a two-phase model to account for empirically derived t 1 values. We present a 
correction of a recent work by Giammanco & Beckman, who use a two-phase model to estimate the ionization rate of H n regions by 
cosmic rays. We show that a very small amount of cold gas is sufficient to account for t 2 values typically inferred for H n regions. 
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1. Introduction 



q , Temperature fluctuations in H n regions are a much-discussed problem. Peimbert d 19671 investigated for the first time the effects 
O ,■ of such fluctuations on temperatures empirically derived from spectroscopic observations and found that they may lead to higher 
electron temperatures being derived from the collisionally excited [Om] nebular-to-auroral forbidden-line ratio, r([Om]), than 
those derived from the Balmer jump of the Hi recombination spectrum, T(Hi). If significant temperature fluctuations exist in Hn 
regions and yet are ignored in the analysis, they may lead to underestimating ionic abundances calculated from collisionally excited 
lines (CELs) (e.g. Esteban et al. 2002 ). The parameter f 2 (Peimbert 1967) was introduced to quantitatively characterize temperature 
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fluctuations. Since then the parameter has been extensively used in nebular studies. Under the conditions that [T(r) - 7o] 2 /T ( 2 <k 1, 



o 

. £h ! where T(r) is the local electron temperature and To the average value weighted by the square of density. The value of t 2 can be 
, determined either by comparing T([0 ra]) and T(H i) or by comparing ionic abundances derived from CELs and from recombination 
' lines (RLs) (see Peimbert et al. 2004, for further details). 

Two-phase models, which approximate a nebula by two components of different physical conditions, represent an over- 
simplified, yet frequently used method of studying nebular physics (e.g. Viegas & Clegg 1994, Zhang et al. 2005). Using a two-phase 
model, Stasinska (2002 ) points out that the classical picture of temperature fluctuations may be misleading under certain conditions, 
and three parameters are needed to characterize temperature inhomogeneities. One of the open questions in the study of H n regions 
is that values of f 2 derived from observations are consistently higher than those predicted by photoionization models (Stasinska 
2000). Recently, Giammanco & Beckman (2005; GB05 thereafter) constructed a two-temperature-phase model capable of explain- 
ing t 2 values deduced for a number of H n regions by means of incorporating a component of cool ionic gas ionized by cosmic rays. 
However, f 2 values deduced from observations cannot be applied directly to two-phase models. This is because empirical values 
of t 2 deduced from observations were calculated from formulae derived by assuming random and small-amplitude temperature 
fluctuations, assumptions that are apparently broken for a two-phase model. 

The purpose of the current work is to quantitatively study the relationship between the t 2 values predicted by two-phase models 
and those measured by observations. We show that, in the case of a very cold ionic gas component embedded in a 'normal' Hn 
region, t 2 deduced from observations using the empirical method may have significantly overestimated the real values. 

2. Analysis 

According to Peimbert (1967), for a given ionic species of number density N,, the thermal structure of an Hn region can be 
characterized by an average temperature To and a mean square temperature fluctuation parameter f 2 , defined as 

f T e N e NjdV 

To = — r (1) 

/ N e N t dV 
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and 



J(T e -T Q ) 2 N e NidV 



T 2 J N e NidV 

respectively. Assuming f 2 <k 1, To and t 2 can be determined from measured r([Ora]) and T(Hi) using relations, 
1/9.13 xlO 4 



r([Om]) = 7b 



and 

rcHi) = r (i 
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(Peimbert 1967; Garnett fl992| l. Esteban et al. ((2002') report t 2 values for a sample of Hn regions. Their values are reproduced in 
the first row of TableQ] Although they were obtained by comparing ionic abundances derived from CELs and from RLs, we assume 
that they are the same as would be deduced from from T([0 in]) and T(H i) using Eqs. (3) and (4). The assumption is supported by 
studies of Torres-Peimbert & Peimbert (2003 ), who find general agreement between t 2 values inferred from the differences between 
T([0 in]) and T(H i) and those inferred from the apparent discrepancies between CEL and ORE abundances. Esteban et al. (2002 ) 
did not provide values of To. For the purpose of comparison, we follow GB05 and assume that To = T([0 m]), as given in the second 
row of Table[T] However, as pointed out previously by Stasinska (2002 ), the empirical method of estimating t 2 using Eqs. (3) and 
(4) may be invalid under certain conditions. In the following, we show that in the two-phase model, values of t 2 that are required 
to account for the measured differences between T([0 in]) and T(H i) (or differences between ORL and CEL abundances) may be 
much lower than that derived from the empirical method. 

In the framework of the two-phase model, which assumes that the electron temperature structure of an H n region consists of a 
hot and a cold phase, we follow GB05 and assume equal densities of the two phases and an ionization fraction of unity for the hot 
gas. Electron temperatures are designated as 7h and T c for the hot and cold phases, respectively. The intensity of an [O in] forbidden 
line transition of wavelength A is given by 



/([Olll]) a ~ J N(0 2+ )N e T; l/2 exp(-AE/kT e )dV, 

where AE is the excitation energy of the upper level. It follows that 



(5) 



/([O 1111)4959. 
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= C x exp 



33000 
T([Om]) 



= Cx 



/V(0 2+ ) h /v e h V h 77 1/2 exp(-29200/T h ) + 2V(0 2+ ) c A£V c :r c ~ 1/2 exp(-29200/T c ) 



7V(0 2+ ) 1 W£v h 7 , ~ 1/2 exp(-62200/T h ) + /y(0 2+ ) c /V e c V<T ( : 1/2 exp(-62200/T c ) 



-1/2, 



(6) 



where C is a constant depending only on atomic data, and super- or subscript 'h' and 'c' refer to quantities of the hot and cold 
phases, respectively. In the cold phase, the ionization fraction of oxygen is expected to be lower than that of hydrogen. As a 
reasonable approximation, we assume that 7Y(0 2+ ) C /7Y(0) C = 0.17Y(H + ) C /7Y(H) C . Our analysis is insensitive to this assumption since 
the cold phase essentially contributes no [O ra] forbidden line fluxes given its very low temperature and the fact that emissivities of 
forbidden lines decline exponentially with decreasing temperature (c.f Eq.01. We thus obtain 



7/([Oiii]) = 33000 hT 



77 1/2 exp(-29200/T h ) + 0.1 0x 2 T; l/2 exp(-29200/T c ) 



7T 1/2 exp(-62200/T h ) + O.\0x 2 T; v1 exp(-62200/T c ) 



(7) 



where, using the notation of GB05, is the mass ratio of the cold to hot gas, and x is the ionization fraction of hydrogen in the cold 
gas. It can be easily seen that for 7h » T c , T([0 in]) ~ T^. 
Similarly, from the flux of the Balmer jump, 



7(Bal,3646)~ j N(H + )N e T; 3/2 dV , 

and the flux of Hj3, 

I(HJ3)~ j N(H + )N e T; 5/6 dV , 



we have 



r(Hi) = 



(T h 3/2 + 0x 2 T; V2 \ 
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Fig. 1. t 2 versus [T([Om]) - 7*(Hi)], deduced from the empirical method [Eqs. © and (|4), dotted lines] and for two-phase model 
[Eqs. 0, (10]), CT2]), and ([T3j, solid lines]. Left panels: Ox 2 = 0.01; Right panels: Ox 2 = 0.1; Upper panels: T h = 10000 K; Lower 
panels: T h = 15000 K. 



or 

2 (T h Y V2 r(Hir 2 / 3 r, 2/3 -i 

0X 2 = \—\ ^ TTT . (11) 

\Tc) i-r(Hi)- 2 / 3 r c 2/3 

Equation (fTOl shows that T(H i) weights towards T c . Following GB05, for a two-phase model, we obtain 
To = \ - ; , (12) 



1 + 6x 2 ' 



t 2 = Ox 2 — - — ^— I , (13) 
6x 2 = (J±Z2^ = , (14) 



T h -T c 



2 



(T - T c ) 2 (To - T c ) 2 



and 



T h = To + —5— . (15) 

Therefore, for the two-phase scenario, To and t 2 should be determined from Eqs. (T7]i, (fTOl . ( fT2l . and ( fT3l instead of from Eqs. <£3j 
and ©. The latter are only valid for random and small amplitude fluctuations. 

In Fig. [ljfor given values of 9x 2 and Th, we compare t 2 as a function of T([0 in]) - 7\H i) derived from the empirical method 
using Eqs. ([3} and ©, and that derived in the scenario of two-phase model using Eqs. Q, ( fTOb , and ( [TBI . The plots show that, 
depending on 9x 2 above a critical value of T([0 in]) - T(H i), the empirical method significantly overestimates f 2 , particularly for 
the case of small Ox 2 . The amount of deviation is insensitive to the value adopted for T^. In addition, we find that as the temperature 
difference between the two phases is larger than a critical value (typically ~ 6000 K), empirical t 2 deduced from observations can 
no longer be used at their face values to constrain two-phase models, a point overlooked by GB05 as discussed in the following 
section. 
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Table 1. Estimated values of f 2 , To, and Ox 2 for T c = 100, 1000, and 4000 K for a sample of Hn regions. 
Top and f 2 are taken from Esteban et al. (f2002). The numbers in parentheses are values derived by GB05 
and included here for comparison. 



Object 


NGC 604 


NGC 5461 


NGC 5471 


NGC 2363 


4 


0.027 


0.041 


0.074 


0.128 


Tq.e 


8150 


8600 


14100 


15700 


r([Orn])(4r 0iE ) 


9052 


9943 


15913 


18528 


r(Hi)(4,7o, E ) 


7783 


8011 


12358 


12344 


(f 2 )ioo 


0.00013 


0.00016 


0.00009 


0.00013 


(f 2 )iooo 


0.00413 


0.00053 


0.00312 


0.00429 


(f 2 )4000 


0.02511 


0.03503 


0.02401 


0.03540 


(To)ioo 


9051 


9941 


15911 


18526 


(To) 1000 


9011 


9884 


15860 


18445 


(Xo)4000 


8518 


9370 


15413 


17719 


(to 2 )ioo 


0.00013(0.028) 


0.00017(0.042) 


0.00010(0.075) 


0.00013(0.130) 


(Ox 2 ) 1000 


0.00522(0.035) 


0.00659(0.052) 


0.00356(0.086) 


0.00480(0.146) 


(0X 2 ) 4OO o 


0.08690(0.10) 


0.10664(0.14) 


0.04378(0.14) 


0.05905(0.23) 



3. New estimates of Gx 2 for GB05 model 

GB05 showed that the ionization of cold neutral gas by cosmic rays may significantly contribute to temperature fluctuations. They 
used a two-phase model to explain f 2 values obtained by Esteban et al. (120021 1 for a number of H 11 regions. However, the high 
temperature difference between the two phases in GB05 model (see their Table 1) suggests that t 2 values obtained by Esteban et al. 
cannot be applied directly to two-phase models. Values of 9x 2 derived by GB05 need to be re-considered. 

We re-estimate Ox 2 values for the sample of Hn regions of Esteban et al. (2002). Following GB05, three values of T c are 
assumed, 100, 1000, and 4000 K. Under these conditions, temperature in the cold gas is too low to collisionally excite the [Oin] 
lines, and consequently Eq. (0 can be simplified to 

7X[OlIl]) = 7h. (16) 

Substituting values of Tg and t 2 given by Esteban et al. (2002) [hereafter referred as Tq E and f 2 , in order to distinguish them from 
those deduced from Eqs. (fT2l and (fT3l ll into Eqs. (0 and (01), we obtain T([0 m])(f 2 ,7o,E) and 7\H i)(t^,Top), which are tabulated in 
Rows 3 and 4 of TableQ] Note that here we take the same assumption by GB05 that T^e = T([0 in]). Then 6x 2 can be determined 
from Eq. (HB, where T(U 1) = T(U i)(f 2 , T ,e), T h = T([0 in])(/|, T ,e)- Although the values of T([0 m])(f|,r ,E) thus obtained are 
slightly higher than the actual values deduced from observations (see Table [T), the resulting 6x 2 are hardly affected. 

In Table Q] we compare our 0x 2 values to those of GB05 (given in parentheses); for low values of 8x 2 differences of up to a 
factor of a hundred are found. It can easily be seen that the discrepancies increase with decreasing temperature of the cold gas, as 
suggested by Fig.[T] For T c = 100 K, our derived values of 6x 2 are very low, suggesting that the values of f 2 reported by Esteban et 
al. (f2002 ) can be explained by the existence of a very small amount of cold gas. 

Table[T]also gives t 2 and To values derived from the Eqs. ( fT2l and ( fT3l . As the Table shows, the real t 2 values are lower when 
T c <sc Tq than those derived from the empirical method. As a result, values of the cosmic ray ionization rate, derived by GB05 
have been grossly overestimated [c.f. their Eq. (17)]. Our conclusion is consistent with the range of values inferred for the Orion 
nebula from Gamma ray observations. 

4. Conclusion 

We have studied the relationship between values of f 2 predicted by a two-phase model and those derived empirically from observa- 
tions (empirical method). Our results show that the existence of extremely cold gas within H n regions may lead to overestimated f 2 
calculated from empirically determined T([0 ill]) and T(H 1). We stress that care should be taken when using the two-phase model 
to study large temperature fluctuations of H 11 regions. In this model, CELs are hardly produced by the cold gas, which on the other 

—3/2 

hand makes a large contribution to the flux at the Balmer jump, due to the T e dependence of 7(Bal, 3646). Accordingly, the 
existence of a very small amount of cold material may lead to a large discrepancy between r([Oin]) and r(Hi). In other words, 
in spite of its small mass, the existence of extremely cold material can reproduce apparently large ? 2 (as derived from the empirical 
method), much larger than the actual value [as defined by Eq. (2)]. 

Finally, we revisited the GB05 study of cosmic ray ionization as a mechanism for creating temperature fluctuations in H 11 
regions. While this provides a potential mechanism for creating cold ionized plasma in H 11 regions, we show that the values of 
£ required to produce the ionization have been overestimated in their treatment, due to the t 2 discrepancy discussed above. Based 
on the formulae presented here, we re-estimated their model parameters. The corrections are apparent, particularly in cases where 
temperature of the cold gas component is low, resulting in lower values of £ that agree better with the estimates for the Orion nebula 
published in the literature. 
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